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Abstract. We have studied the nearby old pulsar PSR B 1929+ 10 and its surrounding interstellar medium utilizing the sub- 
arcsecond angular resolution of the Chandra X-ray Observatory. The Chandra data are found to be fully consistent with the 
results obtained from deep XMM-Newton observations as far as the pulsar is concerned. We confirm the non-thermal emission 
nature of the pulsar's X-radiation. In addition to the X-ray trail already seen in previous observations by the ROSAT and XMM- 
Newton X-ray observatories, we discovered an arc-like nebula surrounding the pulsar. We interpret the feature as a bow-shock 
nebula and discuss its energetics in the context of standard shock theory. 
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1. Introduction 

', In 1990s, a series of powerful X-ray observatories were 
• launched to space. Among many great results on various as- 
' pects in astrophysics, these observatories allowed a first de- 
tailed study of the X-ray emission properties of rotation- 
powered pulsars as a class (cf., Becker & Triimper 1997; 
Becker & Pavlov 2001, for a review). Most olcfl radio pul- 
sars, though, were still too faint for a detailed study by these 
satellites (cf.. Sun et al. 1993; Manning & Willmore 1994; 
Becker & Triimper 1997; Saito 1998). Thanks to the much 
improved sensitivity of the XMM-Newton and Chandra obser- 
[ vatories, a more detailed study of old pulsars became possi- 
ble. Observations of PSR B0950+08; B0823+26; J2043+2740 
(Becker at al. 2004; ZavHn & Pavlov 2004); B0628-28 (Becker 
et al. 2005); B0943+10 (Zhang, Sanwal & Pavlov 2005); 
Bl 133+16 (Kargaltsev, Pavlov & Garmire 2006); B1929+10 
(Becker et al. 2006); and B2224+65 (Hui & Becker 2007a) for 
the first time have allowed us to study the emission proper- 
ties of old pulsars as a class. Surprisingly, the X-ray emission 
from old pulsars seems to be dominated by non-thermal radi- 
ation processes. Thermal components (e.g., to account for the 
emission from hot polar caps) are not required to model the 
X-ray spectra of these pulsars, resulting in mostly upper lim- 
its for any thermal radiation components. Further support for 
an emission scenario dominated by non-thermal radiation is 
found by the observed temporal emission properties. Pulse pro- 



' In standards of high energy astronomy rotation-powered pulsars 
are called young, middle aged, and old if their spin-down age is of the 
order of few times 10^ - 10* yrs, 10' - 10* yrs, and > 10* yrs, respec- 
tively. This classification is dilfuse, though, with a smooth transition 
in between the different groups. 



files, if detected with sufficient photon statistics, show multiple 
components and narrow features. This is indicative of strongly 
beamed emission, which further invalidates the heated polar 
cap scenario as the main source of X-ray emission in old pul- 
sars. The pulsed fractions in old pulsars are in the range of 
~ 30 - 50%. 

In addition to the pulsar emission, which originates within 
the co-rotating magnetosphere, extended trail like X-ray emis- 
sion was observed from PSR B 1929+ 10 (Becker et al. 2006) 
and PSR B2224+65 (Hui & Becker 2007a) on a scale of sev- 
eral arc-minutes. In the case of PSR B 1929+ 10 the X-ray emis- 
sion in the trail is interpreted as synchrotron emission pro- 
duced in the shock between the pulsar wind and the surround- 
ing medium (see the discussion in Becker et al. 2006 for de- 
tails). Owing to the moderate (15" half energy width) spatial 
resolution of XMM-Newton, details of the nebular emission as- 
sociated with PSR B1929+10 might have been remained unre- 
solved. With a ten times improved angular resolution Chandra 
data can thus be essential to further constrain the properties of 
the nebula very near to the pulsar. 

In this paper, we present a detailed analysis of multi-epoch 
Chandra observations of the field around PSR B 1929+ 10. 
Based on its X-ray emission properties, this source is consid- 
ered prototypical of an old pulsar (Becker et al. 2006). With 
a pulse period of P = 226.5 ms and a period derivative of 
P = 1.16 X 10"'^, its characteristic age is determined to be 
~ 3 X 10* years. These spin parameters imply a spin-down lu- 
minosity of E = 3.9 X 10^^ erg s"' and a magnetic field at the 
neutron star magnetic poles of ~ 5 x 10^' G. With a ra- 
dio dispersion measure of 3.178 pc cm"-', the NE2001 Galactic 
free electron density model of Cordes & Lazio (2002) predicts 
a distance of 170 pc. However, the recent astrometric measure- 
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Table 1. Ephemerides of PSR B 1929+ 10 



Right Ascension (J2000) 


19''32'" 13.983' ±0.002" 


Declination (J2000) 


-HlO" 59' 32.41" ±0.07" 


First date for valid parameters (MJD) 


52929 


Last date for valid parameters (MJD) 


53159 


Pulsar rotation period (s) 


0.2265182954 


Period derivative P (10"'^ s s"') 


1.164739 


Characteristic age (10* yrs) 


3.09 


Surface dipole magnetic field (10'^ G) 


0.5129 


Dispersion Measure (pc cm"^) 


3.178 


Distance (pc) 


361!^" 


Spin-down Luminosity (10^^'') ergs s"' 


3.89 



"Adopted from Becker et al. (2006) 



ments by Chatterjee et al. (2004) yielded a precise proper mo- 
tion and parallax determination that translates into an accurate 
distance measurement of d - 36 1;^^" pc and a proper motion 
of y_L = 177+4 g-i jjjg ephemerides of PSR B 1929+ 10, 
which we made use of in this paper, are listed in Table 1 . 

The paper is organized as follows: in Sect. 2 we describe 
the observations and data analysis; and in Sect. 3 we summarize 
the results and provide a discussion in the context of a standard 
shock model for the nebular emission. 

2. Observations and data analysis 

Data analysis is restricted to the energy range 0.5 - 8.0 keV. All 
energy fluxes, however, are computed for the 0.5 - 10 keV band 
for better comparison with the results based on XMM-Newton 
data (Becker et al. 2006). 

The object PSR B1929+10 was observed with Chandra in 
2005 December 04 (Obs ID: 6657) and 2006 May 28 (Obs ID: 
7230) with the Advanced CCD Imaging Spectrometer (ACIS). 
In both observations, PSR B1929+10 was located on the back- 
illuminated (BI) ACIS-S3 chip with an off-axis angle of ~ 0.1 
arcmin. We used standard processed level-2 data. The effective 
exposures are ~ 21 ks and ~ 25 ks for the observations in 2005 
December and 2006 May, respectively. 

2.1. Spatial analysis 

With the objective obtaining better statistics for the analysis, 
we combined both datasets to produce better images. We care- 
fully checked and corrected the aspect offsets for each observa- 
tion prior to the merging. 

The combined X-ray images of the 4x4 arcmin field cen- 
tered on PSR B1929+10 as well as a close-up of the central 
1x1 arcmin regions are shown in Fig. 1. A compact neb- 
ula, which has an arc-like morphology somewhat resembling 
a bow-shock, is clearly detected around PSR B 1929 +10. X-ray 
contours were calculated at the levels of (6 - 26) x 10 counts 
arcsec"^ s"' and overlaid on the image in Fig. lb. Assuming 
this arc-like compact feature is a bow-shock pulsar wind neb- 
ula, it is interesting to estimate the stand-off angle 0^ which is 



a crucial parameter in discussing the shock physics (see also 
Sect. 3). By fitting a 2-D Gaussian model to the raw image of 
the field around PSR B1929+10, we obtained a full width at 
half maximum (FWHM) of 1.13 + 0.02 arcsec (Icr error) . We 
take this as an estimate for the lower bound of 6s- 

To further examine the extent of the compact nebula, we 
computed the brightness profile in the energy band 0.5 - 8 keV 
from the raw image of bin size of 0.5 arcsec. We estimated the 
counts in sixteen consecutive boxes of dimensions 1x15 arc- 
sec, oriented along the nebula. The vignetting coiTected bright- 
ness profile and the orientation of the 16 spatial bins are shown 
in Fig. 2. The arc-like nebula extends up to ~ 10 arcsec from 
the pulsar before fading into the background. Because of the 
wings of the ACIS-S point spread function, the pulsar is ex- 
pected to contribute to its surrounding nebula. To estimate this, 
we compared the observed data with the ACIS-S instrumental 
point spread function. To do so, we simulated a point source 
observation using MARX (ver 4.2. 1). The adopted inputs were 
the pulsar spectrum and the same exposure, roll, and off-axis 
angle as in the actual data. We then computed the brightness 
profile of the simulated data in the same way as aforemen- 
tioned. The estimated point source profile is shown in Fig. 2 
as a dotted line. The brightness profile of the arc-like nebula 
computed from the actual observation is significantly different 
from the expected distribution of a point source. On the other 
hand, there is no significant extended emission detected ahead 
ofPSRB1929+10. 
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Angular distance from PSR B1929+10 (arcsec) 

Fig. 2. Vignetting corrected brightness profile of the field 
around PSR B1929+10. The energy range is 0.5-8 keV. The 
profile expected for a point source is indicated by a dotted line. 
The average background level and its Icr deviation are indi- 
cated by horizontal lines, which were calculated by sampling 
the source-free regions around the pulsar within a 1 x 1 arcmin 
field-of-view. The inset shows the sixteen bins used in comput- 
ing the profile. Each bin is 1" x 15". 

In addition to the arc-like nebula near the pulsar, we ob- 
served a clumpy structure (which is labeled as C in Fig. lb). 
There are only ~ 24 net counts from this clump, which does 
not allow for a more detailed analysis. The signal-to-noise ra- 
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15.0 15.0 19:32:14.0 13,0 12,0 

Right Ascension (J2000) 

Fig. 1. (a) 4x4 arcmin image centered on PSR B1929+10 in the energy band of 0.5 - 8 keV. Top is north and left is east. The 
image is created by merging both Chandra observations. It is binned with a factor of 0.5 arcsec and adaptively smoothed with a 
kernel of cr < 2 arcsec. (b) The close-up view of the central 1x1 arcmin region as illustrated in a. The proper motion corrected 
pulsar position is indicated by the white cross. An arc-like feature associated with the pulsar is observed. The proper motion 
direction of the pulsar (Chatterjee et al. 2004) is illustrated by the arrow. Contour lines at the levels of (6.1 - 26) x 10 counts 
arcsec"- s"' are overlaid. 



tio of this clump is estimated to be ~ 2, and thus is consistent 
with background fluctuations. 

In Fig. la, no prominent structure resembling the X-ray 
trail seen by ROSAT and XMM-Newton opposite to the pul- 
sar's proper motion direction can be identified (cf. Fig. 2 in 
Becker et al. 2006). We have examined the Chandra image fur- 
ther by smoothing the raw data with a kernel of cr < 6 arc- 
sec, which is comparable with the FWHM of XMM-Newton's 
PSF The smoothed image is displayed in Fig. 3. Comparing 
the image with the contours calculated from XMM-Newton 
MOSl/2 data, a faint trail like feature is noticed. However, 
the background contribution of this feature is estimated to be 
~ 60%. The unabsorbed X-ray flux of a 1 arcmin circular 
trail region near the pulsar, as detected by XMM-Newton, is 
/x = 5.3 X lO-i'^ erg/s/cm^ within the 0.5-10 keV band (Becker 
et al. 2006). The low significance of the trail in the two Chandra 
observations, thus, is in agreement with Chandra's lower sen- 
sitivity. We, therefore, will not further consider the X-ray trail 
emission in this paper, and alert the interested reader to the 
XMM-Newton results reported by Becker et al. (2006), for a 
detailed discussion of its emission properties. 

2.2. Spectral analysis 

Although the spectrum from PSR B1929-I-10 has already 
been tightly constrained by the XMM-Newton data (Becker et 
al. 2006) it is important for us to re-examine its spectral prop- 
erties with Chandra since it resolves the pulsar emission from 
the compact surrounding nebular component better 



We extracted the pulsar spectrum in both data sets from 
a circle of 2 arcsec radius (encircled energy ~95%), cen- 
tered on the pulsar, and fitted them simultaneously. The back- 
ground spectrum was extracted from a source-free region 
within a 10 arcsec radius centered at R.A.=19''32™14.107\ 
Dec.=10°59'51.57" (J2000). After background subtraction, 
593 + 24 and 682 ± 26 net counts were available for the spec- 
tral analysis. These values imply the net counting rates of 
(2.84 + 0.1 1) X 10-^ cts/s and (2.77 + 0.1 1) x 10"^ cts/s for the 
observations in 2005 December and 2006 May, respectively. 

We computed the response files with the tools MKRMF and 
MKARF in CIAO 3.4. Utilizing the most updated calibration 
data, CALDB 3.4.1, the generated response files corrected the 
degradation of the quantum efficiency in the ACIS CCD, ac- 
cordingly. Each spectrum was dynamically binned so as to have 
at least 30 counts per bin. To better constrain the spectral prop- 
erties, we fitted simultaneously the spectra obtained from both 
observations. We performed all the spectral fittings in 0.5 - 8 
keV with XSPEC 11.3.2. The parameters of afl fitted model 
spectra are summarized in Table 2. All the quoted errors are 
1 - cr, and were computed for 1 parameter of interest. 

Among the single component models tested, we found that 
a power-law model fits the data best (xl = 0.89 for 36 D.O.F.). 
This model yields a column density of Nh - 2.40+g34 x 10^^ 
cm"^, a photon index of F = 2.91^q }3, and a normalization at 
1 keV of 8.12+i °4 X 10"^ photons keV-' cm^^ g-^ . These best- 
fit values are fully consistent with those obtained by XMM- 
Newton for power-law fits (see Table 3 in Becker et al. 2006). 
The best-fit power-law spectrum and residuals are shown in 
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Fig.3. The 5x5 arcmin region around PSR B1929+10 in the 
energy band of 0.5 - 8 keV. The image is created by merging 
both Chandra observations. It has a binning factor of 2 arcsec 
and is adaptively smoothed with a kernel of cr < 6 arcsec. Top 
is north and left is east. The contours from the XMM-Newton's 
MOS 1/2 0.2 - 10 keV band image at the levels of 0.09 and 0.44 
counts arcsec"^ are overlaid. The proper motion corrected pul- 
sar position is indicated by the white cross. The proper motion 
direction of the pulsar is indicated by the arrow in the image. 

Fig. 4. The column density inferred from both Chandra and 
XMM-Newton spectra is relatively high in comparison with the 
mean value obtained from the Extreme Ultra Violet Explorer 
(EUVE) measurements of the stars in the neighborhood of 
the pulsar (Slowikowska et al. 2005). However, as discussed 
by Slowikowska et al. 2005, the interstellar medium is very 
patchy in the region around PSR B1929-H10 (cf.. Table 3 in 
Slowikowska et al. 2005 for more details). Merely on the basis 
of the comparison with the mean value, the Nh we report and 
in Becker et al. (2006) cannot be invalidated. A more detailed 
mapping is thus needed to better constrain the Nh toward the 
pulsar. 

We have also computed the error contours to demonstrate 
the relative parameter dependences of the photon index vs. the 
hydrogen column density and plotted this in Fig. 5. The unab- 
sorbed flux deduced for the best fit power-law model parame- 
ters is fx - 2.5 X 10"'^ ergs s"' cm"^ within 0.5 - 10 keV. At 
a distance of 361 pc it implies a luminosity of Lx = 3.9 x 10""^ 
ergs s"'. 

As for the XMM-Newton data it is obvious that the single 
power-law model already describes the observed pulsar spec- 
trum very well. Hence, the justification of including an addi- 
tional thermal component is absent. In the Chandra data, we 
found that fitting with a power-law plus blackbody model does 
not yield a reasonable solution if we allow the blackbody radius 



1 2 
channel energy (keV) 



Fig. 4. Energy spectra of PSR B 1929-1- 10 as observed with the 
Chandra ACIS-S3 detector on 04 Dec 2005 and 28 May 2006 
and simultaneously fitted to an absorbed power-law model {up- 
per panel) and contribution to flie;if^ fit statistic (lower panel). 

and the temperature to be free parameters. It resulted in a black- 
body radius of 7?bb = 25.81^'5 g| m. This is much smaller than 
the canonical size of a polar cap, r^^ - RilnRjcPyi^ ~ 300 m, 
which suggests the contribution from the additional blackbody 
component is insignificant. We quantified the statistical signifi- 
cance for adding this extra component to the power-law model 
with the F-test which suggests that inclusion of this thermal 
components is only required at a confidence level of 26%. 

We have also examined the possible contributions from a 
polar cap and the neutron star surface by fixing the blackbody 
radii ?A R - 300 m and R - \0 km, respectively. These fits 
yield temperatures of T ~ 8.7 X 10^ K and T ~ 5.6 x 10^ K. 
The F-test suggests that adding these thermal components to 
the power-law model is only significant at a confidence level 
$ 62%. Due to the better photon-statistics this number was 
even smaller in the XMM-Newton data (Becker et al. 2006). 
This low significance is also reflected by the relative contribu- 
tion of the thermal component in the total energy flux observed 
by Chandra. We estimated the Icr upper limits for the polar cap 
temperature and its flux contribution by adding a blackbody 
component to the power-law model. Fixing the blackbody ra- 
dius at a polar cap size of R - 300 m, a Icr upper limit of 
T < 1 .2 X 10^ K was computed by assuming contribution from 
one polar cap only. This implies that the Icr limit contribution 
from the polar cap in the total energy flux is ~ 12%. 

For the spectral model consisting of two blackbody com- 
ponents, we found that the best-fit model parameters also 
agree well with those inferred from the XMM-Newton spec- 
tra (Becker et al. 2006). Despite the acceptable value of the 
goodness-of-fit, the inferred blackbody radii are too small to 
be considered as a reasonable description. Moreover, the pul- 
sar spectrum obtained by XMM-Newton and fitted together 
with the spectrum obtained by the ROSAT PSPC have already 
shown that this model cannot describe the data beyond ~ 5 keV 
(see Becker et al. 2006). 

Since the arc-like nebula is resolved by Chandra for the 
first time, it is instructive to examine its energy spectrum, de- 
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Table 2. Spectral parameters inferred from fitting the Chandra ACIS-S data obtained from PSR B 1929+ 10 and its associated 
extended feature. 



Model " 




D.O.F. 


Nh 
10^' cm-2 


T" IkT 


Normalization at 1 keV^ 
photons keV"' cm"^ s"' 


Radius '' 
m 


PSR B 1929+ 10 


PL 


0.89 


36 




9 Q1+0.16 
^•^^-0.13 


8.12!i:0^ X 10-5 




BB 


2.00 


36 


<0.04 


0.36!«:«1 




29!^ 


PL+BB 


0.93 


34 


1 40+1"' 


2.48!?:«/0.29!°:g 


4.60:f;J» X 10-5 


25.81!>|^t 


PL+BB 


0.92 


35 


4.4.+0.60 




8.17!;:«xl0-5 


300 


PL+BB 


0.90 


35 


9 79+1.06 
^•'^-0.55 


2 97+0.27 ,QQ5 +0.01 

^•^'-0.17'"-"-^-0.05 


8.67:^;g X 10~5 


10000 


BB+BB 


0.91 


34 


<0.70 


71+0.21/0 27+0.03 

' -0.13' -0.05 




<: ^7+3.43 749 ^9+22.74 
-2.30' -4.47 


Arc-like feature 


PL 


0.88 


6 


2.40 


2 00+0-32 


3.63^^:^5 X 10-6 





" PL = power-law; BB = blackbody 

The entry in this column depends on the model in interest. It is the temperature kT in keV or the photon index T 
" The normalization constant for the power-law model. 

The radius of the blackbody emitting area is calculate for an assumed pulsar distance of 361 pc. 




spite the poor photon statistics. We extracted the nebular spec- 
tra from the regions in accordance with its observed morphol- 
ogy. For the sake of consistency, we adopted the same extrac- 
tion regions in both datasets, which are illustrated in Fig. 6. It 
consists of an annular region centered on the pulsar position 
with inner and outer radii of 2 arcsec and 5 arcsec, a circu- 
lar region with a radius of 2.2 arcsec, as well as a ellipse of 4 
arcsecxS arcsec. 



Fig. 6. 1x1 arcmin raw image (0.5 - 8 keV) centered on 
PSR B1929+10. Data from both Chandra observations were 
used. The extraction regions used for the spectral analysis of 
the arc-like feature are illustrated. 



The background spectra were extracted from each dataset 
within a nearby source-free region of a 10 arcsec radius, 
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channel energy (keV) 

Fig. 7. Energy spectra of the arc-like structure associated with 
PSR B 1929+ 10 as observed by the Chandra ACIS-S3 detector 
on 04 Dec 2005 and 28 May 2006, and simultaneously fitted to 
an absorbed power-law model (upper panel) and contribution 
to the;^'^ fit statistic {lower panel). 



centered at RA =19''32"15.470' and Dec.=+10°59'29.01" 
(J2000). After background subtraction, there were 29 + 5 and 
41+7 net counts extracted from the arc-like feature, imply- 
ing the net counting rates of (1.29 ± 0.26) x 10""^ cts/s and 
( 1 .65 ± 0.28) X 1 0"^ cts/s for the observations in 2005 December 
and 2006 May, respectively. Within the 1 cr errors of these count 
rates, no variability is apparent on the basis of these two obser- 
vations. The response files were computed in the same manner 
as those for the pulsar spectra. Each spectrum was dynamically 
binned so as to have at least 10 counts per bin. In order to obtain 
a better statistic, we analyzed both spectra simultaneously. 

We hypothesized that the nebular emission originates from 
the interaction of the pulsar wind and the ISM. Synchrotron 
radiation from the ultra-relativistic electrons is generally be- 
lieved to be the emission mechanism of the pulsar wind neb- 
ula, which is characterized by a power-law spectrum. We tested 
this hypothesis by fitting an absorbed power-law model to 
the nebular spectra. We fixed the column density at the value 
Nh = 2.4 X 10^' cm"^ inferred from fitting the pulsar spec- 
trum. The power-law model describes the observed spectrum 
reasonably well (xv = 0.88 for 6 D.O.F.). The best-fit power- 
law spectrum and residuals are shown in Fig. 7. This model 
yields a photon index of F 



1 keV of 3.63!°-^^ x lO-*" 



= 2.00;^Q3^ and a normalization at 
photons keV"' cm"^ s"'. The un- 



absorbed flux deduced for the best-fitted model parameters are 
fx = 1.7 X 10"^^ erg s"' cm"^ in the energy range of 0.5 - 10 
keV. The pulsar distance of 361 pc implies a luminosity of 

Lx = 2.7 X lO^'' erg s"'. 

We have checked the robustness of all the spectral parame- 
ters quoted here by incorporating background spectra extracted 
from different source-free regions. We found that within the 
Icr errors the spectral parameters inferred from independent fit- 
tings are aU consistent with each other. 



Table 3. Unabsorbed fluxes, luminosities, and the conversion 
efficiencies of PSR B1929+10 and the nebular components in 
0.5 - 10 keV 



Component 




fx 




a 


Lx/E 




ergs 


-1 -2 

s cm 


ergs 


S-' 




PSR B 1929+ 10 


2.5 


X IQ-" 


3.9 X 


1030 


1.0 X 10-3 


Trail-like feature * 


5.3 


X lO-''* 


8.3 X 




2.1 X 10-" 


Arc-like feature 


1.7 


X lO-''* 


2.7 X 




6.9 X 10-5 



" Luminosities are calculated at the assumed pulsar distance of 361 pc 
* Adopted from Becker et al. (2006) 



3. Discussion 

Complementing the XMM-Newton observation, we have stud- 
ied PSR B 1929+ 10 and its surrounding medium with Chandra, 
which provides data with sub-arcsecond resolution. So far, 
three distinct components have been resolved, namely the pul- 
sar itself, the X-ray trail opposite to the pulsar's proper mo- 
tion direction as well as the arc-like nebula surrounding the 
pulsar. The flux contributions from these components and the 
corresponding X-ray conversion efiiciencies are summarized in 
Table 3. 

The morphology of the arc-like nebula and the orientation 
with respect to the pulsar's proper motion direction suggest 
a bow-shock nature. This is the first weU-defined bow-shock 
morphology observed from an old pulsar. 

A bow-shock nebula as observed in PSR B 1929+ 10 can 
be produced by a pulsar in supersonic motion. In this case, the 
termination shock radius Rs is determined by the balance of the 
ram pressure between the relativistic pulsar wind particles and 
flie ISM at the head of the shock (cf. Hui & Becker 2007b): 



R, ^ iE/2npisMvW'^ ~ 3 X IQi 



16/7l/2„-l/2„-l 



'E--n 



Vp.ioo cm 



(1) 



where Vpjoo is the velocity of the pulsar in units of 100 km 
s"'; £34 is the spin-down luminosity of the pulsar in units of 
10^^ erg s '; and n is the number density of the ISM in units 
of cm"^. A density of 1 cm"^ implies a shock radius of Rg ~ 
10'* cm. The stand-off" angle 6s between the pulsar and edge 
of the shock in the forward direction is related to Rs by 0, - 
Rs cos i/D, where D is the distance to the pulsar and ; is the 
incUnation angle between the proper motion and the plane of 
the sky. Taking the FWHM resulted from Gaussian fitting to 
the radial profile (i.e., 1.13 arcsec) as the lower limit of Gg, i is 
then constrained to be < 50°. 

At a distance of 361 pc, the proper motion velocity and 
the size of the arc-like nebula imply the pulsar took ~ 150 
years to pass through the sky region of the nebula. We es- 
timated the magnetic field in this shocked region by equat- 
ing the passage time to the synchrotron cooling timescale 
Tsyn = enmec/yo-TB^. ^ 10^ (/zv/keV)""-^ B^^^^ yr, where y 
is the Lorentz factor of the wind, taken to be 10^ (cf., Cheng 
et al. 2004), crj is the Thompson cross section, and Bj^q is the 
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magnetic field in the shocked region in units of micro gauss. 
This implies a magnetic field of ~ 75 jiG for the region near to 
the pulsar. For comparison, the typical field strength in the ISM 
is ~ 2 - 6 juG (cf., Beck et al. 2003, and references therein). If 
the magnetic field in the neighborhood of PSR B1929+10 is 
comparable with this estimate, the inferred compression factor 
is ~ 13 - 38 in the bow-shock region, which is rather high 
in comparison with other pulsar wind nebulae (e.g., Hui & 
Becker 2006, 2007b). Polarization measurements of the sky re- 
gion around PSR B 1929-1-10 in the radio band can help to better 
constrain the magnetic field of this system and hence the com- 
pression factor. 

In the context of the standard shock theory (ChevaUer 
2000), the X-ray luminosity and spectral index depend 
on the inequality between the characteristic observed fre- 
quency V^^ and the electron synchrotron cooling frequency: 
Vc = l^nenieC I (T^Tly^B^ (cf.. Chevalier 2000 and references 
therein), which is estimated to be Vc = 1.8 x 10^^ Hz. Since, in 
general v** > Vc, this suggests the X-ray emission is in a fast 
cooling regime. 

We further probed the energy distribution of the syn- 
chrotron radiating electrons in the shock. The emitting elec- 
trons are distributed as A'(7) oc y^P . In a fast cooling regime, 
the theoretical luminosity per unit frequency is given by 
(cf., Cheng et al. 2004): 

Assuming the energy equipartition between the electron and 
proton, we took the fractional energy density of electrons eg to 
be ~ 0.5 and the fractional energy density of the magnetic field 
€b to be ~ 0.01. We then integrated Eq. 2 over the frequency: 
Lf = [""'L^pdv, with hv„,in =0.5 keV and hv„„^ =10 keV. 
Equating L* and the observed luminosity of the arc-like neb- 
ula results in an index of p = 2.35 for the electron spectrum. 
This estimate Ues within the acceptable range of the index p as 
suggested in a standard shock model, i.e., p ~ 2 - 3 (cf., Cheng 
et al. 2004 and references therein). For a fast cooling scenario, 
the photon index is related to p as T'^ = ip + 2)/2. The value 
p - 2.35 implies a photon index off* = 2.18. This theoretical 
value lies within a Icr uncertainty of the observed value. 

The arc-shaped morphology of the nebula observed in 
PSR B1929+10 is similar to that of Geminga (Caraveo et 
al. 2003). However, there is a main difference between these 
two cases. While the arc-like nebula of Geminga is rather sym- 
metric with respect to the direction of proper motion, asymme- 
try is indicated in the case of PSR B 1929-1- 10 (see Fig. lb). Its 
asymmetric shape might be a result of the anisotropic pulsar 
wind and/or inhomogeneities in the surrounding ISM (see Hui 
& Becker 2006). 

For the pulsar PSR B1929-I-10, we found that the spec- 
tral properties inferred from our analysis are in good agree- 
ment with the results obtained by XMM-Newton, which did 
not allow the pulsar emission to be separated from the com- 
pact nebular component. The consistency is expected even in 
the presence of the diffuse compact nebulae as it contributes 
only ~ 7% to the observed energy flux (see Table 3). Thus, 



the non-thermal emission scenario of PSR B1929-f-10 is con- 
firmed. Hui & Becker (2007a) had argued that it is possible to 
sustain particle acceleration regions in PSR B 1929-H lO's outer- 
magnetosphere (so-called outer-gap) if the inclination of the 
magnetic axis with respect to the rotational axis is taken into 
account. This inference is supported by the abihty of the outer- 
gap model in reproducing the observed X-ray pulse profile and 
its phase shift relative to the radio pulse (see Fig. 16 Becker et 
al. 2006). However, this model has difliculty in explaining the 
observed spectral properties of PSR B1929-I-10. 

In the outer-gap model, the non-thermal X-rays result from 
the back-flowing charge particles from the outer gap (Cheng 
& Zhang 1999). When the primary electrons/positrons leave 
the outer-gap, they will emit curvature photons that are sub- 
sequently converted into secondary pairs in the presence of 
the strong magnetic field. Synchrotron photons will then be 
emitted by these secondary electrons/positrons. If these pho- 
tons are energetic enough, they will further be converted into 
pairs that again lose their energy via synchrotron radiation. 
Therefore, an electromagnetic cascade is developed. Based on 
this model, Cheng & Zhang (1999) argued that the X-ray pho- 
ton index resulting from such cascades should be < 2. This is 
obviously not in agreement with the observed photon index for 
PSR B1929-H10, which is as steep as 2.91![J J^. Observations 
of five other old pulsars: B1133H-16 (Kargaltsev,Pavlov & 
Garmire 2006); B0943-I-10 (Zhang, Sanwal & Pavlov 2005); 
B0628-28 (Becker et al. 2005); B0823-H26 (Becker et al. 2004); 
and J2043H-2740 (Becker et al. 2004), also found the photon in- 
dices steeper than 2. This gives the outer-gap emission model a 
challenge. Re-examination of the model is thus required. 

After this work was submitted to A&A for publication, 
we became aware that a paper on the similar subject as pre- 
sented here was submitted to ApJ by Misanovic, Pavlov, & 
Garmire (2007). The results reported by these authors agree 
well with the results presented here. Misanovic et al., though, 
report to have detected a thermal component in the spectrum of 
PSR B1929-I-10. While their parameters inferred from the sin- 
gle power-law fittings are consistent with ours, they suggested 
that a model of power-law plus blackbody model provides a 
statistically better description. This obviously contradicts with 
the conclusion presented here and in Becker et al. (2006). From 
a closer inspection of their spectral fits we conclude that their 
apparent detection is probably biased by a too small signal- 
to-noise ratio per spectral bin. While we binned the Chandra 
spectrum so as to have at least 30 counts per bin Misanovic et 
al. (2007) grouped the data so as to have only about 15 counts 
per spectral bin. For the spectral fitting of XMM-Newton data, 
the parameters obtained by Misanovic et al. (2007) are for 
190 degree of freedom whereas Becker et al. (2006) had only 
121 degrees of freedom for their analysis. Becker et al. (2006) 
thus used a much higher photon statistics per spectral bin. The 
higher signal-to-noise ratio per spectral bin yielded a more 
stringent discrimination among competing spectral models. For 
this reason, the spectral results presented here as well as in 
Becker et al. (2006) provide more stringent constraints on the 
pulsar's X-ray emission properties than those in Misanovic et 
al. (2007). 
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